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Abstract 

Introduction: The recent availability of efficacious prevention interventions among stable couples offers new opportunities for 
reducing HIV incidence in sub-Saharan Africa. Understanding the dynamics of HIV incidence among stable couples is critical to 
inform HIV prevention strategy across sub-Saharan Africa. 

Methods: We quantified the sources of HIV incidence arising among stable couples in sub-Saharan Africa using a cohort-type 
mathematical model parameterized by nationally representative data. Uncertainty and sensitivity analyses were incorporated. 
Results: HIV incidence arising among stable concordant HIV-negative couples contribute each year, on average, 29.4% of total 
HIV incidence; of those, 22.5% (range: ll.l%-39.8%) are infections acquired by one of the partners from sources external to the 
couple, less than 1% are infections acquired by both partners from external sources within a year and 6.8% (range: 3.6%-11.6%) 
are transmissions to the uninfected partner in the couple in less than a year after the other partner acquired the infection from 
an external source. The mean contribution of stable HIV sero-discordant couples to total HIV incidence is 30.4%, with most of 
those, 29.7% (range: 9.1%-47.9%), being due to HIV transmissions from the infected to the uninfected partner within the 
couple. The remaining incidence, 40.2% (range: 23.7%-64.6%), occurs among persons not in stable couples. 
Conclusions: Close to two-thirds of total HIV incidence in sub-Saharan Africa occur among stable couples; however, only half of 
this incidence is attributed to HIV transmissions from the infected to the uninfected partner in the couple. The remaining 
incidence is acquired through extra-partner sex. Substantial reductions in HIV incidence can be achieved only through a 
prevention approach that targets all modes of HIV exposure among stable couples and among individuals not in stable couples. 

Keywords: stable couples; sources of infection; HIV incidence; Sub-Saharan Africa; demographic and health surveys; 
mathematical model. 
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Introduction 

The majority of the adult population in sub-Saharan Africa 
(SSA) live in marital or co-habiting couples [1]. A considerable 
fraction of these stable couples (SCs) are HIV sero-discordant 
(that is, one partner testing HIV sero-positive while the other 
testing HIV sero-negative) [1-5]. Recent randomized clinical 
trials have demonstrated the potential for averting much 
of the HIV heterosexual transmission among stable HIV sero- 
discordant couples (SDCs) using highly efficacious HIV preven- 
tion interventions such as antiretroviral therapy (ART) [6-9]. 
This progress in HIV prevention research revealed new 
horizons for HIV policy and programming in SSA and placed 
SDCs on the list of priorities for HIV prevention efforts [10-12]. 
Yet, reaching a consensus on an effective HIV prevention 
strategy in SSA that factors SDCs entails a comprehensive 
understanding of the sources of HIV incidence arising in the 
population, particularly among SCs. 



We recently quantified the contribution of HIV incidence 
within SDCs to total HIV population-level incidence (SDC, nt ) 
across SSA [13]. We also estimated the fraction of HIV 
infections arising among SDCs that are due to sources 
external to the couple [14]. We further quantified the risk 
of HIV transmission from the infected to the uninfected 
partner in the SDC [15]. In this article, we complement 
our work by extending our analysis to the sources of HIV 
incidence across the different types of SCs for 24 countries 
in SSA. Specifically, we assess the contribution of HIV 
incidence among stable concordant HIV-negative couples 
(SCNCs) relative to total HIV population-level incidence 
where: a) one of the partners acquires the infection from a 
source external to the couple (SCNC extxl ), b) each of 
the partners acquires the infection from a source external 
to the couple (SCNC extX 2) and c) one partner acquires 
the infection from a source external to the couple and 
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then transmits the infection to the uninfected partner in 
the couple shortly after acquiring the infection (SCNC e)rt+int ). 
We also estimate the contribution of HIV incidence among 
SDCs where the uninfected partner acquires the infec- 
tion from a source external to the couple (SDC e)tt ). We 
further provide an estimate for the incidence contribution 
of individuals not in a SC where HIV infection is acquired 
through sexual contacts outside of the context of marital or 
cohabiting partnerships (NSC). Finally, we update the results 
of our earlier work on SDCs [13,14] through the use of 
recently available survey data, in addition to including 
countries that were not covered earlier, and the use of 
a more elaborate mathematical modeling framework for 
generating these estimates. 

Methods 

We constructed a cohort-type mathematical model that esti- 
mates the contribution of HIV incidence stratified by couple 



status and source of infection to total HIV incidence in the 
population over the course of a year. The model uses the 
calculated risks of HIV transmission and acquisition among 
SCs to derive the annual number of new HIV infections arising 
among SCs and the relative contribution of each source of 
incidence to the total HIV incidence in the population. 

Model structure and measures' definition 

Our model sets out to ask the following question: if there is a 
national screening survey at a point in time (at Time 0) that 
identifies individuals engaged in SCs and others not engaged 
in SCs, what proportion of these individuals would acquire 
HIV over the following year (at Time 1), that is within a year 
of observation, and how? Figure 1 shows the possible 
changes to the HIV sero-status of the three types of SCs 
that are identified in the baseline cross-sectional survey at 
Time 0 (SCNCs, SDCs or concordant HIV-positive couples), 
and changes to the HIV sero-status of individuals not in SCs, 
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Individuals not in a 
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Individuals not in a stable couple acquiring HIV infection through 
sexual contacts outside the context of marital or cohabiting 
partnerships. 





Figure 1. Model conceptualization for HIV incidence in the population classified based on the sero-status of stable couples and source of 
infection. The table shows the possible outcome scenarios and the associated mathematical expressions for the different HIV incidence 
measures. The green circle indicates an HIV sero-negative individual, while the red circle indicates an HIV sero-positive individual. 

*Parameters include k: the probability of an HIV sero-negative partner in a stable couple (SC) to acquire the infection from a source external 
to the couple over the course of one year; N couples : the number of SCs identified in the baseline screening cross-sectional survey at Time 0; 
Pscnc- the prevalence of stable concordant HIV-negative couples among all couples; P SDC : the prevalence of stable HIV discordant couples 
among all couples; t 6mths : the probability that the index partner who acquired the infection from an external source will transmit the infection to 
the uninfected partner during the six months following the acquisition of HIV; t lyear : the probability that the index partner in a stable HIV 
discordant couple will transmit the infection to the uninfected partner during the time between the two cross-sectional surveys at Time 0 and 
Time 1; N rep age : the size of the population in reproductive age; ) ' ln _ coup i es : the fraction of the population in reproductive age engaged in SCs; 
P: HIV prevalence in the population; ip: the HIV population-level incidence rate. 
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one year later at Time 1. The mathematical expressions used 
to assess the annual number of new HIV infections arising 
among these population sub-groups are also summarized in 
this figure and discussed further in the Supplementary file. 

In summary, we drew a map of HIV incidence stratified 
by couple status and source of infection by measuring, for 
each country, the contribution of six different types of HIV 
incidence to total HIV incidence arising in the population 
over the course of a given year of observation (Figure 1). 
The latter is defined as the number of new HIV infections 
arising within a year among susceptible individuals in reproduc- 
tive age. This measure defines the denominator in all six 
contribution measures, and is calculated using HIV population- 
level incidence rate and the estimated number of uninfected 
individuals in reproductive age for each country (Table 1 and 
Supplementary file). 

Contributions among SCNCs: 

1) SCNC extxl is the contribution of new HIV infections 
arising among SCNCs where only one of the partners 
acquires the infection from a source external to the 
couple, relative to total HIV incidence in the population. 
The number of infections here is calculated using the 
probability of acquiring the infection from a source 
external to the couple by one of the partners and the 
number of SCNCs in the population (Figure 1 and 
Supplementary file). 

2) SCNC extx2 is the contribution of new HIV infections 
arising among SCNCs where both partners acquire the 
infection from a source external to the couple within 
the same year, relative to total HIV incidence in the 
population in that year. The number of infections here is 
calculated using the probability of acquiring the infec- 
tion from a source external to the couple by both 
partners and the number of SCNCs in the population 
(Figure 1 and Supplementary file). 

3) SCNC ext+ , nt is the contribution of new HIV infections 
arising among SCNCs where the index partner acquires 
the infection from a source external to the couple, 
on average half the way through the year of observa- 
tion, and transmits HIV to the other partner during the 
following half of the year, relative to total HIV incidence 
in the population. The number of infections here is 
calculated using the probability of acquiring the infec- 
tion from a source external to the couple by one of the 
partners, the likelihood of transmitting the infection to 
the uninfected partner within six months following HIV 
acquisition and the number of SCNCs in the population 
(Figure 1 and Supplementary file). 

Contributions among SDCs: 

1) SDC, nt is the contribution of new HIV infections arising 
among SDCs where the HIV infected partner transmits 
the infection to the uninfected partner, relative to 
total HIV incidence in the population. The number of 
infections here is calculated using HIV transmission 
probability per partnership over the course of a year of 



observation and the number of SDCs in the population 
(Figure 1 and Supplementary file). 
2) SDC e)rt is the contribution of new HIV infections arising 
among SDCs where the uninfected partner acquires 
HIV from a source external to the couple, relative to 
total HIV incidence in the population. The number 
of infections here is calculated using the probability 
of acquiring the infection from a source external to the 
couple and the number of SDCs in the population 
(Figure 1 and Supplementary file). 

Contribution among individuals not in SCs: 

1) NSC is the contribution of new HIV infections arising 
among susceptible individuals in reproductive age who 
are not part of a SC, relative to total HIV incidence 
in the population. The number of infections here is 
calculated using HIV population-level incidence rate and 
the number of susceptible individuals who are in 
reproductive age but not part of a SC (Figure 1 and 
Supplementary file). 

Model parameterization 

Multiple data sources were used to obtain the model 
parameters with the primary source being the Demographic 
and Health Surveys (DHS), which are standardized nationally 
representative household-based surveys [16]. We analyzed 
the most recent DHS data for 24 countries in SSA where an 
HIV serological biomarker survey has been conducted [16]. 
These data were complemented by population size informa- 
tion from the United Nations World Population Prospects 
Database [17] to calculate country-specific demographic, 
behavioural and epidemiological indicators (Table 1). 

As per DHS, a SC is defined as a man and a woman living 
in a consensual union within a household at the time of 
the DHS cross-sectional survey [18]. Accordingly, polygamous 
arrangements may contribute multiple SCs. SCs where one or 
both partners did not test for HIV were excluded from 
our analysis. Missing HIV information among all SCs ranged 
from 0.5% to 27.3% (mean of 10.8%) across countries. DHS 
guidelines were followed in applying to our calculations the 
sampling weights retrieved from the DHS datasets [18,19]. 
Further details related to the management of DHS datasets 
can be found in the Supplementary file. 

Following the methodology applied in our earlier work 
[13], we calculated the total HIV incidence in the population 
using HIV population-level incidence rates estimated using 
the Joint United Nations Programme on HIV/AIDS (UNAIDS) 
SPECTRUM model [20]. When these were not available, 
HIV population-level incidence rates were derived from the 
DHS HIV prevalence measures using the expression HIV 
population-level incidence rate = ^population prevalence [21] Further 

~ r Duration ot infection L 1 

details can be found in the Supplementary file. 

To estimate HIV incidence arising among susceptible 
individuals that are not in SCs, we assumed that the risk of 
HIV acquisition among these individuals is equal to the aver- 
age HIV population-level incidence rate among all susceptible 
individuals in the population (Figure 1 and Supplementary 
file) [14]. Country-specific values of the HIV population-level 
incidence rates can be found in Table 1. 
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Table 1. Key demographic and HIV-related indicators across the 24 sub-Saharan African countries included in our analysis 
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Countries are shown in order of increasing HIV prevalence in the population. 

Pop: population; Rep: reproductive; Prev: prevalence; SDC: stable HIV discordant couples; MC: male circumcision; Inc: incidence that is the number of new HIV infections per year. 
*Data on male circumcision were not collected during the 2010-2011 round of the DHS for Senegal. The rates used are drawn from a previous Senegal DHS survey conducted in 2005; ¥ HIV population-level 
incidence rate estimated by SPECTRUM or derived using DHS HIV prevalence per 100 person-years; 3f: mean probability of acquiring HIV from sources external to the couple per 100 person-years derived 
by performing 10,000 runs of model fits. 
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The risk of HIV transmission within SDCs was calculated 
using the best available empirical evidence for HIV transmis- 
sion probability per heterosexual coital act (p) as measured 
in the Rakai Study [22] and in the Partners in Prevention 
HSV/HIV Transmission Study (Partners in Prevention Study) 
(Table 2) [23,24]. To adjust for male circumcision (MC) and 
self-reported condom use at last sex among SCs, we varied 
the efficiency of HIV transmission among SDCs reporting 



one or both interventions by applying to p multiplicative 
factors adjusting for the efficacy of MC and condom use in 
preventing HIV transmission. The effect of MC was stratified 
based on whether the index partner in the couple is a male 
or a female. We incorporated the coverage of condom use 
through a fraction of the acts that are protected by the 
partial efficacy of this intervention. The mathematical expres- 
sions and further details on the calculation of the risk of HIV 



Table 2. Model assumptions in terms of key parameter values related to HIV transmission and acquisition in sub-Saharan Africa 



Assumptions 



Parameter values 



Source 



Probability of acquiring HIV from sources external to the 

couple per 100 person-years (X) 
Probability of acquiring HIV by an individual not in a stable 

couple per 100 person-years (<p) 
HIV transmission probability per coital act (p) 

Acute infection (p ocute ) 

Latent infection [pi atent ) 

Average (p) using the Rakai and the Partners in 
Prevention Studies 
Average (p) using the Rakai Study 
Average (p) using the Partners in Prevention Study 
Frequency of coital acts per month (n) 
Demographic and epidemiological measures 

Number of individuals in reproductive age in the 
population (N repage ) 

Number of stable sexual couples identified in baseline 
screening survey (N couples ) 

Fraction of the population in reproductive age that are 

engaged in stable couples (f lncouples ) 

HIV prevalence in the population (P) 

Prevalence of stable concordant HIV-negative couples 

(Pscnc) 

Prevalence of stable HIV discordant couples (P S dc) 
Fraction of females (index partners) among those initially 
concordant HIV-negative couples (/,„d exS cNc) 
Fraction of females (index partners) among 
HIV discordant couples {f in d e xsac) 
Fraction of circumcised males among concordant HIV- 
negative couples \f mcp0 p) 

Fraction of circumcised males among HIV discordant 
couples where the female is HIV infected {f mcin dex) 
Fraction of coital acts protected by condom use (/ m „d om ) 

Efficacy of condoms in preventing HIV transmission per 
condom-protected coital act (E condom ) 

Efficacy of male circumcision in preventing HIV acquisition 
among males per coital act (f mc ) 

Duration 

Between each round of the cross-sectional survey 

{^follow-up) 

Acute infection (t acute ) 

Latent infection spent by index partner between two 
subsequent cross-sectional surveys (T fate „ t ) 



Derived from model fits 

HIV-population-level incidence rate estimated by 
SPECTRUM or derived using DHS HIV prevalence 

0.036 
0.0008 
0.00115 

0.0012 
0.0011 
8.3 acts per month 

Table 1 

Table 1 

Table 1 

Table 1 
Table 1 

Table 1 
50% 

Table 1 

Equal to fraction of circumcised males in the 
population (Table 1) 
Table 1 

Table 1 
80% 

58% 



1 year 

49 days 
134 days 



Derived 
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transmission in presence of these interventions can be found 
in the Supplementary file. 

The likelihood of an HIV sero-negative partner in a SC 
acquiring the infection from a source external to the SC 
(defined here as A.) was determined by the condition that 
the total HIV incidence in the population, as estimated 
from pooling together all incidence measures among SCs and 
individuals not in SCs, must be equal to the total HIV 
incidence as estimated independently using the SPECTRUM 
model or derived using the DHS HIV prevalence (Table 1). 
All incidence measures were then recalculated using this 
fitted value of A, and estimates for the contributions of new 
HIV infections among SCs to total HIV incidence in the 
population were derived (Figure 1 and Supplementary file). 

Uncertainty and sensitivity analyses 

For each country, uncertainty analyses were performed for the 
estimates of the contribution of each type of HIV incidence 
by implementing 10,000 runs of the model using Monte 
Carlo sampling from triangular probability distributions for 
the ranges of demographic, biological and epidemiological 
parameters (Table SI in the Supplementary file). Parameter 
ranges were primarily determined by the 95% confidence 
intervals (CIs) around the country-specific DHS measures. 
In the event where CIs were not available to provide a range, 
plausibility ranges informed by the range of available data 
in the literature or general consensus in the field were used 
as parameter ranges. In each run, model fits were conducted 



Individuals not in 
stable couples 
(NSC) 
40.2% 



and an estimate for the probability of acquiring HIV from 
sources external to the couple (A) was derived. Country- 
specific distributions for the estimated contributions of 
SCs to total HIV incidence by couple status and source of 
infection were then generated and used to calculate the 
mean and associated 95% CIs of these estimates (Figures 3 
and 4). 

Sensitivity analyses were conducted, using Kenya as an 
example, to assess the sensitivity of the calculated contribu- 
tion measures to variations in level of condom use at last sex 
among SCs, MC coverage in the population and the fraction 
of the population in reproductive age that are engaged in SCs 
(Figure S2 in the Supplementary file). 

Results 

The key demographic and HIV-related indicators for the 
24 countries in SSA can be found in Table 1. The estimated 
numbers of new HIV infections per year among adults for 
each country are also included in Table 1. Figures 2, 3 and 4 
show the contributions of the six measures of HIV incidence 
stratified by couple status and source of infection. Over the 
course of a year of observation, stable concordant HIV- 
negative couples (that is SCNCs) contribute on average 29.4% 
of the total HIV population-level incidence across these 
countries in SSA (Figures 2 and 3). Most new HIV infections 
among these couples are due to one partner acquiring the 
infection from a source external to the couple (SCNC e)rtxl ), 
with an average of 22.5% (range: ll.l%-39.8%) across all 



External HIV 
acquisition by both 
partners in a 
concordant 
negative couple 
within the same 
year (SCNCext*2) 
Wk. 0.07% 




External HIV 
acquisition by one 
of the partners in a 
concordant 
negative couple 
(SCNCext*1) 
22.5% 



HIV transmission shortly after the 
external acquisition of HIV by the 
index partner in a concordant 
negative couple (SCNCext+int) 6.8% 



External HIV 
acquisition by the 
uninfected partner 
in a discordant 
couple (SDCext) 
0.7% 




HIV sero-discordant couples (30.4%) 



Figure 2. The average contributions to the total number of new HIV incident infections in a year in the population stratified by couples' 
sero-status and source of HIV infection for 24 countries in sub-Saharan Africa. The average for each mode of exposure represents an average 
over the country-specific mean contribution measures (fraction of new HIV infections relative to total HIV incidence in the population in a 
given year). For each country, the mean contribution of each source of exposure to total HIV incidence was calculated based on 10,000 runs 
of the model using Monte Carlo sampling from triangular probability distributions for the specified ranges of model parameters. 
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♦ Mean contribution of stable concordant HIV-negative couples where one of the partners acquires the 
infection from a source external to the couple (SCNCext'1 ) 
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♦ Mean contribution of stable concordant HIV-negative couples where each of the partners acquire 

the infection from a source external to the couple (SCNCext*2) 
Mean SCNCext*2 across all countries 
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♦ Mean contribution of stable concordant HIV-negative couples where one partner acquires HIV from 
a source external to the couple and then transmits the infection to the uninfected partner in the 
couple shortly after acquiring the infection (SCNCext+int) 
Mean SCNCext+int across all countries 
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Figure 3. Mean and 95% confidence interval of the contributions of HIV incidence among stable concordant HIV-negative couples to total 
HIV incidence in the population in 24 countries in sub-Saharan Africa. The figure shows the contribution of HIV incidence among stable 
concordant HIV-negative couple where: (A) one partner acquires the infection from a source external to the couple, (B) each of the partners 
acquire the infection from a source external to the couple and (C) one partner acquires the infection from a source external to the couple and 
then transmits it to the uninfected partner in the couple. Estimates were calculated based on 10,000 runs of the model for each country using 
Monte Carlo sampling from triangular probability distributions for the specified ranges of uncertainty of the model parameters. Countries are 
shown in order of increasing HIV prevalence. The horizontal line in the different panels represents the average for the contribution measure 
in question across all countries. 
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♦ Mean contribution of identifiable HIV incidence among stable HIV discordant couples due to HIV 

transmission from the infected to the uninfected partner in the couple (SDCint) 
Mean SDCint across all countries 
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contacts outside the context of marital or cohabiting partnerships (NSC) 
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Figure 4. Mean and 95% confidence interval of the contributions of: (A) identifiable HIV incidence among stable HIV discordant couples due 
to HIV transmission from the infected to the uninfected partner in the couple, (B) HIV incidence among stable HIV discordant couples due to 
acquiring the infection from a source external to the couple and (C) HIV incidence among individuals not in stable couples. These measures, 
for 24 countries in sub-Saharan Africa, are relative to total HIV incidence in the population in each country. Estimates were calculated based 
on 10,000 runs of the model for each country using Monte Carlo sampling from triangular probability distributions for the specified ranges of 
uncertainty of the model parameters. Countries are shown in order of increasing HIV prevalence. The horizontal line in the different panels 
represents the average for the contribution measure in question across all countries. 
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Figure 5. Correlation with HIV prevalence of the mean contribution of: (A) stable concordant HIV-negative couples where one partner 
acquires the infection from a source external to the couple (SCNC extvl ), (B) stable concordant HIV-negative couples where each of the 
partners acquire the infection from a source external to the couple (SCNC ext ., 2 ), (C) stable concordant HIV-negative couples where one 
partner acquires the infection from a source external to the couple and then transmits it to the uninfected partner in the couple 
(SCNC ext+/nt ), (D) identifiable HIV incidence among stable HIV discordant couples due to HIV transmission from the infected to the uninfected 
partner in the couple (SDC,„ t ), (E) HIV incidence among stable HIV discordant couples due to acquiring the infection from a source external to 
the couple (SDC ex ,) and (F) HIV incidence among individuals not in a stable couple (NSC). Values for the Pearson correlation coefficients (r) 
and their associated p-values are incorporated. The analysis discounts the uncertainty in these measures (arising from uncertainty analyses). 
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countries (Figure 3A). The average contribution of both 
partners acquiring the infection from an external source 
(SCNC extx2 ), within the same year, is less than 1% (range: 
0.004%-0.4%) across the countries (Figure 3B). The contribu- 
tion attributed to the index partner transmitting the infection 
to the uninfected partner shortly after HIV acquisition from 
an external source (SCNC ext+ , nt ), that is within the same year 
of observation, is on average 6.8% (range: 3.6%-11.6%) 
across countries (Figure 3C). 

New HIV infections occurring among SDCs identified as 
sero-discordant at the onset of the year of observation 
contribute on average 30.4% of the annual total HIV 
population-level incidence across the countries (Figure 2). 
Most of these HIV sero-conversions are due to acquiring the 
infection from the HIV sero-positive partner in the couple 
(SDC, nt ; Figure 4A) with an average contribution of 29.7% 
(range: 9.1%-47.9%). HIV acquisitions from sources external 
to the couple among SDCs (SDC ext ) contribute minimally 
to total HIV incidence with an average of 0.7% (range: 
0.1%-3.1%) across these countries (Figure 4B). A substantial 
proportion of HIV incidence occurs among individuals not 
in SCs (NSC), with a mean of 40.2% (range: 23.7%-64.6%) 
across countries (Figures 2 and 4C). 

Our findings suggest a strong dependence on HIV pre- 
valence for SCNC extx2 (Pearson correlation coefficient (r): 
0.81, p<0.001, Figure 5B) and SDC ext (r. 0.93, p<0.001, 
Figure 5E), where higher contributions are observed in high 
HIV prevalence countries. NSC is moderately correlated with 
HIV prevalence (r: 0.53, p-value: 0.008, Figure 5F). Mean- 
while, there is no evident correlation with HIV prevalence 
for SCNC extxl (r: -0.35, p-value: 0.093, Figure 5A), 
SCNCgxt+^t (r: —0.23, p-value: 0.283, Figure 5C) and SDC, nt 
(r. -0.33, p-value: 0.115, Figure 5D). The 95% CIs around 
our estimates for the contributions, generated using the 
uncertainty analyses, confirmed our finding of roughly equal 
contribution of SCNCs, SDCs and individuals not in SCs to 
total HIV incidence in the population across SSA (Figures 3 
and 4). 

Discussion 

In this article, we complement our earlier work [13-15] 
by conducting a comprehensive quantitative mapping of 
the contributions of all sources of HIV exposure to total HIV 
incidence in the population stratified by couple status. 
For completeness, we also include updates to few measures 
reported earlier [13,14] by analyzing more recent, or first- 
time available, DHS data for a number of countries and 
further improving on the precision of some estimates using a 
more complex mathematical model. Our findings show that 
close to two-thirds of total HIV incidence every year in SSA 
occur within the context of marriage or cohabitation, but 
only half of these infections are actually attributed to HIV 
transmissions from the infected to the uninfected partner 
within a couple. The rest of the infections are acquired 
through extra-partner sexual encounters to the couple. HIV 
incidence in SSA appears to be roughly equally distributed 
among concordant HIV-negative couples, HIV sero-discordant 
couples and individuals not in SCs, with no dominant mode 
of exposure through which individuals acquire HIV (Figure 2). 



Among concordant HIV-negative couples, most HIV inci- 
dence is due to one of the partners in the couple acquiring 
HIV from an external source, while the incidence arising from 
external HIV acquisitions by both partners within the same 
year is minimal. HIV transmission to the uninfected partner 
shortly after the external acquisition of HIV by the other 
partner in the SCNC, is rather limited at less than 10% of HIV 
incidence. Meanwhile, HIV sero-conversions among SDCs 
due to the transmission of the infection from the HIV sero- 
positive to the HIV sero-negative partner in the couple 
contribute the majority of HIV incidence among SDCs. 

These findings are in close agreement with those mea- 
sured in the Rakai Study cohort where the contributions of 
SCNCs, SDCs and individuals not in SCs were 38.5%, 24.0% 
and 37.6%, respectively (after excluding partnerships with 
incomplete HIV sero-status information) [33]. Our findings 
are also in agreement with a recent modeling study that has 
assessed the incidence contributions among SCs [34]. This 
study used a mechanistic model that tracks HIV incidence 
in the population in 18 countries in SSA, as opposed to our 
programming-oriented approach that is based on a func- 
tional definition of incidence among couples in a framework 
of repeated cross-sectional surveys (Figure 1). Though there 
are technical differences in the classification of couples' 
status and the modes of exposure between the approach 
of Bellan et al. [34] and that of our study, the findings of both 
studies converge on the conclusion that HIV incidence in the 
population is distributed roughly in equal proportions among 
external infections to couples, within couples and among 
persons not in couples. 

Our updated estimates for HIV incidence among SDCs 
reaffirm our earlier published work where HIV transmissions 
from the infected to the uninfected partner in an SDC 
contributed about a third of new HIV infections arising in the 
population [13], and the vast majority of HIV incidence 
arising among SDCs [14]. Indeed, our updated findings for the 
risk of acquiring the infection externally among SDCs confirm 
our earlier results as upper-bound estimates for this risk [14]. 

While SCNCs, SDCs and individuals not in SCs appear 
to contribute equally to total HIV incidence in the popula- 
tion across SSA, heterogeneity can still be observed across 
countries (Figure SI in the Supplementary file). One determi- 
nant of the observed differences is the variability in HIV 
prevalence across Africa (Table 1). This is manifested in the 
strong association of SCNC ext x 2 (Figures 3B and 5B) and SDC e)(t 
(Figures 4B and 5E) with HIV population prevalence. Such 
results are expected as both measures depend proportionally 
on X, the likelihood of acquiring the infection from a source 
external to the couple (mathematical expressions in Figure 1), 
which is higher in high HIV prevalence countries compared 
to low HIV prevalence countries. The dependence of X on 
HIV population incidence rate (that is indirectly on HIV 
population prevalence) also explains the absence of a correla- 
tion with HIV prevalence for SCNC ext x 1 (Figures 1, 3A and 5A) 
and SCNC e)rt+ , nt (Figures 1, 3C and 5C). Meanwhile, the 
moderate correlation with HIV prevalence of the fraction of 
the population that are engaged in SCs, which appears to 
be lower in high HIV prevalence countries, is reflected as 
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a moderate correlation between NSC and HIV prevalence 
(Figures 1, 4C and 5F). 

Other factors could also affect the variability of the con- 
tribution measures across SSA. These include the level of 
condom use and coverage of MC in addition to the rates of 
engagement in SCs across countries (Supplementary file). 
Condom use and MC within SDCs reduce HIV transmission 
within SDCs, therefore leading to lower SDC int and implicitly 
higher contributions from HIV infections acquired through 
sources external to the couple to explain the total observed 
HIV incidence (Figures S2A and S2B in the Supplementary 
file). Conversely, higher levels of engagement in SCs are asso- 
ciated with an increase in the contribution of HIV incidence 
occurring among SCs (Figure S2C in the Supplementary file); 
so are the increased prevalence of HIV discordancy and 
HIV-negative concordancy (Figure 1). The impact of varia- 
tions in these factors on contribution measures is magnified 
whenever more than one factor is at play. For example, in 
Swaziland and Lesotho, the world's largest HIV epidemic 
centres, the higher levels of reported condom use among 
SCs (Table 1), the lower rates of engagement in SCs (Table 1) 
and the low prevalence of discordancy among partnerships 
affected by HIV [1] have yielded lower SDC, nt contribution 
(Figure 4A), and lower than expected SCNC extxl (Figure 3A) 
and SCNC ext+ , nt (Figure 3C) contributions but higher NSC 
contribution (Figure 4C). Despite the observed heterogeneity, 
the overall picture appears to suggest an equal contribu- 
tion of SCNCs, SDCs and individuals not in SCs to total HIV 
incidence in the population. 

One of the highlights of our study is quantifying the 
contribution of SCNC e)rt+ „, t . Although, in principle, these 
infections can be considered as occurring among SDCs 
instead of SCNCs, the individuals who are newly infected 
from a source outside the couple are unlikely to be detected 
within a year by typical counselling and testing programs. 
We label this type of HIV incidence as non-identifiable 
HIV incidence among SDCs to distinguish it from that 
where the infected partner transmits HIV to the uninfected 
partner among couples that have been already identified 
as SDCs in a cross-sectional survey before HIV transmission 
to the partner (that is SDC, nt ). The latter is labelled as 
identifiable HIV incidence among SDCs and is discussed at 
length in an earlier publication [13]. Our estimates assume 
annual cross-sectional surveys, but these may not be feasible 
to implement in resource-limited areas. A lower frequency 
of cross-sectional surveys would increase the fraction of non- 
identifiable HIV incidence among SDCs at the expense of the 
fraction of identifiable HIV incidence among SDCs, making it 
harder to implement effective SDCs-targeted interventions. 

It is worth mentioning that our estimates for the con- 
tribution measures factor in implicitly the role of polygamous 
partnerships in driving HIV incidence among SCs. The DHS- 
based samples of SCs include each union between a man and 
a woman at the time of the cross-sectional survey irrespec- 
tive of the number of concurrent unions. 

Our results show that over one-third of HIV population- 
level incidence occurs among individuals not in SCs (Figures 2 
and 4C). This finding is probably not surprising since close to 
half of HIV incidence in SSA occurs among youths [35-37]. 



HIV infections arising among young individuals are probably 
less likely to have been acquired in the context of spousal 
partnerships and are more likely to have been acquired 
by other modes of exposure such as through casual or 
commercial sex encounters. 

Although another third of new HIV infections in the popu- 
lation arises from extra-marital partnerships, determining 
the context in which these infections are acquired, whether 
through contact with commercial sex networks, other forms 
of high risk behaviour, or through casual sex, requires an 
extended mathematical model that includes multiple risk 
groups and heterogeneity in risk behaviour. Such an analysis 
is beyond the scope of this publication whose purpose is to 
quantify the contribution of the different sources of HIV 
incidence occurring among SCs, rather than investigating the 
actual drivers of the HIV epidemic across different settings. 

The nature of our results where HIV incidence is dis- 
tributed over different sources of exposure limits our ability 
to develop specific HIV policy and programming recommen- 
dations for SSA and suggests the need to simultaneously 
address the different modes of HIV exposure in the popula- 
tion. A treatment as prevention (TasP) approach [38], and the 
new World Health Organization treatment recommendations 
of increasing ART coverage to include HIV infected individuals 
with CD4 count <500 [39], should achieve substantial 
reductions in HIV transmission across all sources of exposure 
to HIV infection. However, the specific impact of such inter- 
ventions on HIV incidence through each source of exposure 
is not yet known. ART coverage may vary from one sub- 
population to another due to variable access to voluntary 
counselling and testing services. Maintaining HIV infected 
individuals in the treatment cascade may also vary across 
sub-populations [40]. The impact of the interventions on the 
onward transmission of HIV in the population is also not 
yet entirely clear. With the increasing availability of data on 
ART coverage, more complex mathematical models can be 
designed to elucidate the impact of different HIV interven- 
tions on HIV incidence through each of these modes of 
exposure. 

Our estimates are affected by the representativeness and 
precision of available data. The recency of the epidemiolo- 
gical and behavioural parameters used in our analyses was 
determined by the availability of HIV biomarker information 
in the DHS surveys. However, DHS surveys are increasingly 
being implemented in SSA and only few surveys of those 
used in our analysis were older than five years [16]. For some 
countries, data were only available in small sample sizes 
resulting in wide CIs around some of the measures. Our 
estimates could also have been affected by inherent biases 
in the DHS data such as the variability in response rate to HIV 
testing [41], and selection bias in restricting our analysis to 
couples with complete HIV sero-status information. Despite 
these limitations, DHS surveys are the only standardized 
nationally representative surveys conducted consistently 
across the countries of SSA [16]. 

There are heterogeneities in risk behaviour and not all of the 
population in reproductive age is necessarily sexually active. 
Nevertheless, the measures that we used to parameterize 
our model are averages for the specific sub-populations of 
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interest regardless of the heterogeneities within. Heterogene- 
ities can affect the estimates for a specific sub-group of a sub- 
population of interest, but they may not affect the average 
for the whole sub-population. This is especially true because 
the focus of our study is on the short-term HIV incidence 
within a year of observation where non-linear effects that are 
sensitive to variability are not at play. Long-term projections, 
however, including HIV onward transmission, could potentially 
be affected by variability in risk behaviour. 

Epidemic type or phase can affect the distribution of HIV 
infection in a population, and potentially our results. None- 
theless, our model uses empirical cross-sectional data of 
HIV distribution and discordancy for a specific year, and 
epidemic type and phase implicitly affects and drives these 
data. Therefore effects due to epidemic type or phase should 
be implicitly accounted for in our calculations. 

We parameterized our model using the HIV transmission 
probability per coital act (p) as measured in the Rakai Study 
[22] and the Partners in Prevention Study [23,24] which, 
to date, are the best available empirical evidence for this 
measure. Yet, different biological or behavioural factors 
may affect p across settings. We have already accounted in 
our calculations for the coverage of MC and the uptake of 
condom use among SCs in each country (Table 1). Still, other 
factors may affect p but are difficult to adjust for due to data 
limitations such as the presence of sexually transmitted 
infections other than HIV [42,43], other co-infections that 
increase HIV viral load [44,45], viral factors [46-48], differ- 
ences in viral sub-types that may lead to longer period 
of elevated viral concentration in the early stages of HIV 
infection [48] and host genetics and immunology [46]. Our 
recent analysis of the risk of HIV transmission within SDCs 
in SSA suggests substantial variability across countries [15]. 
Incorporation of such variability can affect the reported 
estimates by increasing the contribution of within-couple 
transmission, wherever the risk of HIV transmission is higher, 
and reducing the within couple HIV transmission, wherever 
the risk of HIV transmission is lower. Similarly, if the country- 
specific HIV population incidence rates were actually lower 
than those estimated by the SPECTRUM model or derived 
using the DHS HIV prevalence, this would increase the con- 
tribution of within couple HIV transmission at the expense of 
other contributions. 

We also parameterized our model using p during acute 
infection as derived by Pinkerton using the Rakai cohort data 
[22,25]. This value is similar to those reported by Hollings- 
worth et al. [49], Powers et al. [50] and Boily et al. [51]. Given 
the small contribution of SCHC ext+jnt (Figures 2 and 3C), it is 
not likely that the uncertainty in p during acute infection [52] 
will affect our findings. 

In the absence of country-specific empirical measures, we 
assumed that the risk of HIV acquisition among individuals 
not in SCs is equal to the HIV population-level incidence rate. 
This assumption is reasonable on balance of our knowledge 
of HIV epidemiology among persons in SCs and not in SCs, 
and the fact that individuals not in SCs constitute a large 
fraction of the population in reproductive age, nearly half of 
it, across countries. Empirical data, such as those of the Rakai 
Study [33], as well as HIV incidence age-distribution patterns 



in SSA that are skewed towards young age [53], support also 
the plausibility of this assumption. 

We structured our model to derive the probability of 
acquiring HIV from sources external to the couple (A) and 
hence, increase the precision of our contribution measures. 
Our estimates for the average X as calculated using 10,000 
runs of model fits were reasonable in terms of their scale, 
since, as expected, they were close yet smaller than the HIV 
population-level incidence rate estimated through SPECTRUM 
or derived from DHS. 

The structure of our model does not consider variations in 
partnership duration among SCs. However, the long durations 
of stable partnerships [15], and the large rates of engagement 
in SCs across SSA (Table 1), suggest that partnership forma- 
tion and dissolution within the course of a year of observa- 
tion would not be substantial enough to affect our findings. 

Finally, we have relied on a volume of data sources derived 
using different methodologies to draw a comprehensive 
assessment for the contribution of SCs to the HIV epidemic 
across the African continent. This may have potentially led to 
inconsistencies that could impact our predictions. Though the 
above mentioned limitations may have affected the precision 
of our quantitative results, our uncertainty analyses suggest 
that they are not likely to affect our findings that no single 
mode of exposure among SCs contributes the majority of HIV 
incidence that is occurring in SSA. 

In conclusion, we presented a comprehensive mapping 
of the contribution of SCs to HIV incidence in SSA. Our 
estimates were based on a mathematical model para- 
meterized by state-of-the-art empirical and nationally- 
representative population-based data. We conclude that no 
single mode of exposure among SCs contributes to the 
majority of HIV incidence. Accordingly, a multi-focus HIV 
prevention strategy that optimizes the use of available 
prevention interventions by targeting the different modes 
of exposure to HIV in the population, each according to 
its contribution weight, is needed to address the HIV 
epidemic in SSA. Further translational research and inter- 
vention impact assessments are needed to delineate more 
specific recommendations for HIV policy and programming. 
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